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Abstract 

Synchronization in an array of mutually coupled systems with a finite time-delay in coupling 
is studied using Josephson junction as a model system. The sum of the transverse Lyapunov 
exponents is evaluated as a function of the parameters by linearizing the equation about the 
synchronization manifold. The dependence of synchronization on damping parameter, coupling 
constant and time-delay is studied numerically. The change in the dynamics of the system due 
to time-delay and phase difference between the applied fields is studied. The case where a small 
frequency detuning between the applied fields is also discussed. 
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Time-delayed systems are interesting because the dimension of chaotic dy- 
namics can be made arbitrarily large by increasing the time-delay and hence 
find applications in secure communications. Most of the chaos based communi- 
cation techniques use synchronization in unidirectional drive response system. 
A limitation which arises in this case is that messages can be sent only in one 
direction. Thus for a two way transmission of signals, a bidirectional coupling is 
required. In this work we deal with an array of three mutually coupled Joseph- 
son junctions with a finite time-delay in coupling and study the system dynamics 
in the presence of an external driving field. The effect of phase difference and a 
small frequency mismatch between the driving fields on a chaotically synchro- 
nized system with time delay is then studied. 



I. INTRODUCTION 



The equation for Josephson junction (JJ) when treated within the Stewart-McCumber 
model is identical to the equation for a driven damped pendulum which has been studied 
theoretically for several routes to chaos[l[ Q]. Since the work done by Belykh, Pedersen and 
Soerensen on chaotic dynamics in Josephson junction (JJ) jjj, 4], a great deal of investigations 
have been carried out as it is an interesting candidate in the field of non linear dynamics. A 
detailed review of chaos in J J can be found in [5|. Time-delay has been studied experimen- 
tally in J J transmission lines represented by a sequence of overdamped underbiased JJs 
The synchronous variation of all the bias currents causes a change in the time required for 
the pulses to pass through the line, and hence possible to provide a required time delay. 

A generalized stability theory for synchronized motion of coupled oscillator systems was 
developed by Fujisaka and Yamada jjlsj]. Since Pecora and Carroll 9] reported the synchro- 
nization of chaotic systems, different types of synchronization such as complete, generalized 
and phase synchronization of chaotic oscillators have been described theoretically and ob- 
served experimentally in different physical situations 10|, llll] • Intermittent synchronization 



has been reported in bidirectionally coupled chaotic Josephson junction [12J. Chaos syn- 
chronization in coupled J J using active control techniques has also been studied 13j. We 
have studied the effect of phase difference between applied fields on suppression of chaos and 
synchronization in parallely coupled JJs [14]. Time-delay is ubiquitous in physical, chem- 
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ical and biological systems due to finite transmission times, switching speeds and memory 
effects. Anticipa tory , la g, p rojective and phase synchronizations have been reported in time- 



delayed systems 
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17l | . However the addition of a time-delay considerably complicates 
stability analysis. The condition for stability of synchronization for unidirectionally cou- 
pled piece- wise linear time-delayed system can be obtained using the Lyapunov-Krasovskii 
method 
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191 ]. Intermittent anticipatory, intermittent lag and complete synchronization 



were found to exist at the same time in unidirectionally coupled nonlinear time-delayed sys- 



tem having two different time-delays 20j. Important functional information about the brain 



states could be obtained by modelling it with equivalent scale-free small-world networks 



2l| . A configuration for such a model was given by considering three bidirectionally cou- 
pled oscillators in a line. The outer elements were seen to synchronize isochronally while the 
middle element lagged behind. It was also verified for a neuron model with three elements 
that the outer elements gets synchronized while remaining lag synchronized with the middle 



one 



221 ] . The phenomena of synchronization of the outer lasers in a system of three lasers 
has been explained using the ideas of generalized synchronization [23] . Isochronal synchrony 
in coupled semiconductor and fiber ring laser models with mutual delay-coupling was also 



studied 



24] 



In this work, we chose three mutually coupled Josephson junctions with time-delay in 
coupling as a model for our studies. To our knowledge, mutually coupled systems described 
by a second order differential equation with time-delay has not been studied earlier for phase 
effects on synchronization. In section II, the model is introduced and section III deals with 
the stability analysis. The effect of phase difference and frequency detuning among the 
applied fields on the system is discussed in section IV. In section V, we discuss the results 
and present a conclusion. 

II. THE MODEL 

The normalized equation of a single Josephson junction as represented by the resistively 
and capacitively shunted (RCSJ) model is {25J 

<j> + f3(j) + sin0 = i dc + i cos(Oi), (1) 
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where is order parameter phase difference, (3 — -| 



2ei c C 



is the damping parameter, i dc 



is the applied dc biasing and i cos(f2t) is the applied rf-field. The dynamical equations for 
three Josephson junctions with a time-delay in coupling can be given as 

01 + + sin 0i = i dc + i cos(fit) - a s 0i - 4>2 (t - t) 

02 + (3<i>2 + sin 02 = « s 0i(t - r) + 3 (t - r) - 20 2 

03 + /^03 + sin 03 = i dc + i cos[(Q + AQ)t + 8} - a s 3 - (j> 2 {t - r) , (2) 

where r is the time delay applied and a s is the coupling term which takes into account of 
the contribution from the delay circuit. AQ is the frequency detuning and 9 is the phase 
difference among the applied fields. Phase synchronization has been studied in a similar 
system of chaotic rotators without time-delay [26j. In the present work we have considered 
the systems to be identical and Eq. (T2]) may be written in equivalent form as 



01 = 01 

ipi = -/3-01 - sin 0i + i dc + i cos(fit) - a s [V>i - if) 2 (t - r)] 

02 = ^2 (3) 
ip 2 = -/#02 - sin 02 + a s [-01 (t - r) + 3 (t - r) - 20 2 ] 

03 = 03 

^3 = -/303 - sin 03 + i dc + i cos[((Q + AQ)t + 6)} - a s [0 3 - 2 (* - 7")] • 

For the case Af2 = ^ = 0, it can be seen from Eq. [2]that the subsystem consisting of the outer 
JJs possess symmetry with respect to interchange of variables and hence may possess iden- 
tical solutions. This type of situations where identical solutions exists for coupled systems is 
known as complete synchronization. Chaotic systems also exhibit complete synchronization 
for some values of parameters. An array of JJs with no delay in coupling is found to be 
chaotically synchronized for the parameter values j3 = 0.3, iq = 1.2, Q = 0.6, i dc = 0.3 with 



0.37 



271 ] . We have selected these values of parameters for numerical studies unless 



specified otherwise. In the presence of a time-delay in coupling, it is found that the dynam- 
ics changes between periodic and chaotic motions. However synchronization is found to be 
unaffected by the time delay. 
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III. STABILITY ANALYSIS 



In order to check for the stability of synchronization and its dependence on various 
parameters, we need to know the transverse Lyapunov exponents (TLE) and its dependence 
on the parameters. The necessary and sufficient condition for stability of the synchronous 
solution is that all the transverse Lyapunov exponents (TLE) calculated with respect to the 
perturbation out of the synchronization manifold should be negative. However calculation 
of Lyapunov exponents gets complicated when time-delays are involved. By linearizing the 
equation for the outer JJs about the synchronization manifold, we can arrive at a necessary 
condition for synchronization, i.e., the sum of the TLE should be negative. If the sum 
of the TLEs is negative, it implies that the phase space is shrinking. Let <f>(t) and ip(t) 
represent the synchronous solution and we define new variables A<f>i(t) = <f>i(t) — <p{t) and 
Ai[)i(t) = ipi{t) — ip(t) with % = 1,3. Here A0j(£) and Aipi{t) are the perturbations of 
the outer oscillators from the synchronization manifold. Linearizing Eqj3] transverse to the 



synchronization manifold , we get after dropping the subscripts 23] 




1 \ / A0 

-P-a.) \A^j 

We have approximated sin(A0) rs A<p as A0 is small. Due to the delay in coupling, 
perturbations will not affect the coefficient matrix until t — 1\ > 2r. The Wronskian of the 
linearized system can be related to the trace of the matrix by Abel's formula 



Wit) 



exp 



jf (-a, - P)(ttj . 



A(p Aip 
A<p Aip 

The Wronskian gives the phase space dynamics of the system. Taking the natural log of the 
Wronskian we get 

\n[W(t)] = In \A<pAip - A^A0| = - / {a s + (3)dt. (4) 

Jti 

This is a monotonically decreasing function of t which means that the phase space volume 
of the system perturbed from the synchronization manifold contracts as a function of time. 
The sum of the transverse Lyapunov exponents is given as 

2 

\\ 

t— >oo t 



2 1 

> A,- = lim -ln\A^Atp - At/A0|, (5) 



The sum of the transverse Lyapunov exponents can be now approximated using Eq. H] as 



Ai + A 2 + (6) 

which is negative indicating that the phase space of the coupled system shrinks to a trajectory 
representing the synchronous state. The sum of the Lyapunov exponents depends on the 
coupling term and damping parameter. Even though the sum of the conditional Lyapunov 
exponents is negative, if one of the exponents is positive, the solution will blow up along the 
unstable direction. This kind of a situation, where isochronally synchronized solution get s 
unstable even though the sum of Lyapunov exponents is negative was addressed in ref. 24]. 

The quality of synchronization is usually quantified using the correlation coefficient (CC) 
given by 

cc = dMt) - (Mm w£> - mm) (7) 

\/(\Mt)-(Mt))\ 2 )\/(Mt)-m))\ 2 ) 

with i = 2,3. The value of correlation coefficient lies between —1 < CC < 1 with large 
value of | CC | meaning better synchrony. The cross correlation of the dynamics is shown 
in FigJjTa) for the JJ system, from which we observe that for very low values of coupling 
constant there is loss of synchrony. In Fig{T](b) the transverse Lyapunov exponents are 
plotted. It can be seen that the sum of the TLE will always be negative as given by EqEJWe 
observe that corresponding to the values where cross correlation is lost, there is a small 
positive value for one of the TLE. Figj2] (a) shows the cross correlation between the outer 
junctions for various values of damping parameter. It can be observed from Figj2] (b) that 
the middle junction remains uncorrelated with the outer ones for most of the values of 
damping parameter. Thus the middle junction mediates synchronization between the outer 
junction while remaining unsynchronized from both the outer junctions. 

In the presence of a time-delay r, the dynamics of the system changes considerably. It is 
observed that for some values of time delay (for r nearly equal to 0.35 onwards) the system 
exhibits periodic synchronized motion. Fig. E^a) is the bifurcation plot for various time- 
delays which shows periodic and chaotic behavior for different delay times. With Fig. H](a), 
we show that the system remains synchronized for most of the values of time-delay. The 
TLEs plotted in Fig. H](b) shows that one of the TLE has positive value for regions where 
CC is lost. 
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IV. EFFECT OF PHASE DIFFERENCE AND FREQUENCY DETUNING 



A phase difference between the applied fields was found to supress chaos effectively in an 
array of coupled Josephson junctions. In this section we study the effect of phase difference 
in bidirectionally coupled time-delay systems. From Eq. [2J we can see that in the presence 
of an applied phase difference, the equation for the outer junctions is no longer identical. 
In terms of the difference variable ip^[ 3 = ipi — 1P3 we can write the equation for the outer 
junctions from Eq. [3] as 

i>i,3 = -Wi,3 ~ 2 cos(^) sin(^) + i' sm(nt + 9/2) - a s ^ 3 (8) 

where 0^3 = <pi — 03, <pf 3 = 0i + 03 and i' Q = 2i sin(#/2) and it can be seen that for a finite 
value of phase difference the outer junctions cannot get synchronized. 

The phase difference which usually desynchronizes the system may be used to suppress 
chaos in time-delayed systems. From Figs. [3] and HI it can be seen that for a time delay 
of r = 0.25 the system is chaotic and synchronized. Fig. [5] shows that an applied phase 
difference of 9 = 0.1 5tt desynchronizes the system with r = 0. But when a time-delay along 
with a phase difference is applied the system exhibits periodic motion. From the time series 
plotted in Fig. [6J it can be observed that the application of a phase difference 9 = 0.15-7T 
along with the time-delay, the system goes to a periodic state. 

To complete the discussion, we consider the case of frequency detuning as it is difficult to 
apply a second frequency which is exactly same as the first. Small deviations in the applied 
frequencies are inevitable. So we consider the system with a small frequency detuning and 
the initial phase difference is taken to be zero. In order to investigate the difference in the 
dynamics brought about by frequency detuning, we first check the case where AQ = in 
Eq|2j i.e, the case without any detuning. Fig. [7] shows the temporal dynamics of the system 
considered with no detuning. Fig. [7(a) shows the time series plot for ipi with no detuning. 
Fig. [7(b) shows the time series plot for the difference between the variables of the outer and 
the middle junction while Fig. [TJ(c) is time series plot for the outer junctions. The outer 
junctions are perfectly synchronized while remaining uncorrelated with the inner junction. 

When a frequency detuning is applied, it can be observed that the outer junctions show 
synchronization in regular intervals. Fig. [8](a) is the time series plot for the variable for the 
outer junction with AQ = 0.004 and it can be seen that a periodic modulation has appeared 
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due to detuning. Here we observe that though the qualitative behavior is repeated, the 
trajectory of the system cannot be completely repeated due to the chaotic segments in the 
evolution process. This type of motion is referred to as breathers. Fig. WJd) shows the 
time series plot for the difference between the variables of the outer and middle junctions 
and it can be seen that they remain uncorrelated. Fig. [S^c) shows the time series plot for 
the difference between the variables of the outer junctions. It can be seen that they get 
synchronized periodically with a period of T = 2ir/AQ. The qualitative behavior of the 
system is not affected by time delays. 



V. RESULT AND DISCUSSION 



In this work we deal with bidirectionally coupled time-delayed systems and study the 
effect of delay and phase difference between the applied fields on synchronization. The 
sum of the Lyapunov exponents transverse to the synchronization manifold is evaluated and 
it is found to be negative indicating the shrinking of phase space of the coupled system 
to a trajectory representing the synchronous state. However cross correlation coefficient 
reveals positive TLE for some values of coupling constant and damping parameter where 
synchrony is lost, though the sum would be still negative. Transverse Lyapunov exponents 
are evaluated numerically and are found to be in good agreement with the analytic results. 
By varying the time-delay, we have analyzed the behavior of the system and it is observed 
that for small values of time delays, periodic motion occur and as the delay time is increased 
chaotic motion reappears. However the system remained synchronized for most of the values 
of time-delay and hence may find applications in secure communication. The study on a 
configuration of three oscillators in a line is of importance because of the current recognition 
that it has resemblance with neuron models. Phase difference between the applied fields 
together with the time delays may be effectively used to suppress chaos. The practical 
situation where a small frequency detuning will be present between the applied fields is also 
studied. As the frequency detuning is like a time dependent phase difference, periodic and 
chaotic motions are repeated with a period of T = 2tt/AQ. Time-delay does not change the 
effect of frequency detuning. Experimentally it is possible to apply a required delay in JJ 
6|. Hence by suitably adjusting the time-delay, chaos may be controlled in JJ devices like 
voltage standards, detectors, SQUIDS etc where chaotic instabilities are least desired. 
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FIG. 1: (a)Crosscorrelation (CC) between the two outer JJs for different values of coupling strength 
with all other values (3 = 0.3, io = 1.2, = 0.6, id c = 0.3 (b) Transverse Lyapunov exponents plotted 
for the same parameter values. r=0.1 in both cases. 



10 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



0.5 



H o 



-0.5 



innnirf rTT i 



(b) 



0.1 0.2 0.3 0.4 



FIG. 2: (a) CC between the two outer JJs for different values of damping parameter with a s = 0.37 
(b)CC between the outer and middle junction. 
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FIG. 3: (a)Bifurcation diagram for ipi against r (b) Bifurcation diagram for ip2 against r. Parameter 
values are (3 = 0.3, i = 1.2, $7 = 0.6, i dc = 0.3, a s = 0.37 
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FIG. 4: (a)Cross correlation for various values of r with the other parameter values as (3 = 0.3, iq = 
1.2, ft = 0.6, i dc = 0.3, a s = 0.37 (b) TLEs plotted for different values of r 




FIG. 5: The outer and inner junctions are desynchronized on the application of a phase difference 
6 = 0.157T with t = 0. 
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FIG. 8: The temporal dynamics of the variable with Af2 = 0.004 for a = 0.37 
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